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We present a measurement of CP asymmetry using a time-dependent Dalitz plot analysis of 
B° — > tv + tv~tv° decays based on a 414 fb -1 data sample containing 449 x 10 6 B_B pairs. The data 
was collected on the T(4S) resonance with the Belle detector at the KEKB asymmetric energy 
e + e~ collider. Combining our analysis with information on charged B decay modes, we perform a 
full Dalitz and isospin analysis and obtain a constraint on the CKM angle <j>2, 68° < 4>2 < 95° as 
the 68.3% confidence interval for the (j>2 solution consistent with the standard model (SM). A large 
SM-disfavored region also remains. 

PACS numbers: 13.25.Hw, 11.30.Er, 12.15.Hh 



In the standard model (SM), CP violation arises from 
an irreducible phase in the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix [l|, Snyder and Quinn pointed out 
that a timc-dcpcndcnt Dalitz plot analysis (TDPA) of 
the decay B° — > pit — > 7r + 7r~7r° 0] offers a unique way to 
determine the angle <p2 [H in the CKM unitarity triangle 
without discrete ambiguities, which cannot be obtained 
from analyses of other modes sensitive to cf>2 such as B — > 
7T7T or pp Q. The TDPA uses isospin and takes into 
account a possible contamination from 6 — s- c? penguin 
transitions. In addition, using measurements of B + — > 



p + n° and p°7r + provides further improvement of the (j>2 
determination [y, |7|. 

In this Letter, we present the result of a TDPA in 
B° — > it + tt~tt decays and a constraint on We use 
a 414 fb -1 data sample that contains 449 x 10 6 BB pairs 
collected on the T(4S l )resonance. The data were taken 
at the KEKB collider [1] using the Belle detector 0. 

In the decay chain T(45) -> B°B° -> (7r+7r-7r 6 )/ta g , 
where /t ag is a final state that distinguishes B° and B°, 
the time- and Dalitz plot-dependent differential decay 
rate is 



oT 

dAt ds+ds- 



\At\/ 



Tb0 {(\A 3w \ 2 + |A3^ 2 )-<ztag-(|A 37r | 2 - |^| 2 ) cos(Am d At)+ 9tag .2Im 



3in(Am d At)} 



Here, A 3t: is the Lorentz-invariant amplitude of the 
B°(B°) — > 7r + 7r~7T° decay, gtag is the 6-flavor charge 
(<?tag = +1 (—1) when / ta g is a B° (B°) flavor eigen- 
state), and At is the decay time difference of the two B 
mesons (t^ — ttag)- The parameters p and q define the 



mass eigenstates of neutral B mesons as pB° ± qB°, with 
an average lifetime r^o and mass difference Amj. The 
Dalitz plot variables s+, s_, and so are defined as 

«+ = (p++Po) 2 , s- = (P-+Po) 2 , s = (p++p-) 2 , (2) 
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where p + , p_, and po are the four- momenta of the 7r + , 
ir~ , and 7r°, respectively, in the decay of B° — > ir + ir~ir°. 

The amplitudes A 37T have the following Dalitz plot 
dependences 



A 3n (s + ,s.) = Yl Us+,s-)A K , (3) 

«=(+,-,0) 

£ / K ( S+)S _)^, (4) 



-^4 3t (s+, s_) 



K = ( + ,-,0) 



where A K (A K ) are complex amplitudes corresponding to 
B°(B°) -> p+tt-^-tt+^Ott for re = +,-,0. Here we 
neglect possible contributions to the B° — > 7r + 7r~7r° de- 
cay other than that of B° — > (/Cot) — > 7r + 7r~7r° and take 
account of them as systematic uncertainties. The func- 
tions f K incorporate the kinematic and dynamical prop- 
erties of B° — ► (ptt)° decays and can be written as 

{ 7Us + ,s-) = T$ =1 { F } :( Sk ) (k = +,-,0), (5) 

where Tj =1 and C F^(s K ) correspond to the helicity dis- 



tribution and the lineshape of p K , respectively. The line- 
shape is parameterized with Brcit-Wigncr functions cor- 
responding to the jo(770) and its radial excitations: 



< ^( s )=BW p{770) 



{ /3« BW p( i45o) 



( 7«BW 



p(1700) 



(6) 

where the amplitudes /3 K and ( j K (denoting the relative 
sizes of two resonances) are complex numbers. We use the 
Gounaris-Sakurai (GS) model 10( for the Breit-Wigner 
shape and the world average for the mass and width 
of each resonance. Though /? B and '7^ can be different 
for each of six decay modes of B°(B°) — > (pir)° in general, 
we assume no such variation, i.e., f3 K = f3 and '7^ = 7, 
in our nominal fit, and address possible deviations from 
this assumption in the systematic error. This assumption 
leads to the relation / K (s+,s_) = / K (s+,s_). 

With this relation and Eqs. ((3|) and ((H), the coefficients 
of Eq. Q are 



E i/« 

*6{+ ,-,0} 



^ E 

K<cr£{+,- 



(Rc[/ K /;]C/± CT Rc - Im[/ K /;]c/±> Im ) 



,0} 



(7) 



Im [ -Al„A 3 „ 



E l/-l 2/ -+ E (Rc[/ K /;]/^ + im[/ K /;]/^) 



kG{+,-,0} 



k<ct£{ + ,-,0} 



(8) 



with 



U± = 



\A K \ 2 ± \A K \ 2 



lm[A R A K *] , 

Rc(Im) [A K A a * ± A K A a *] , 
Re(Im) p K A CT *-(+)3 CT A K * 



(9) 
(10) 

(11) 
(12) 



In 

rr±,Re(Im) 

^ KG 

rRe(Im) 

KG 

The 27 coefficients (J9J) (Tj"2"j) are the parameters deter- 
mined by the fit [T3|. The parameters ©- ((101) and (TTT1) 
(I12p arc called non-interfering and interfering parame- 
ters, respectively. This parameterization allows us to 
describe the differential decay width as a linear combi- 
nation of independent functions, whose coefficients are 
fit parameters in a well behaved fit. We fix the overall 
normalization by requiring Uf = 1. Thus, 26 of the 27 
coefficients are free parameters in the fit. 

In contrast to a quasi-two-body CP violation analysis, 
a TDPA includes measurements of interfering parame- 
ters, which are measurements of CP- violating asymme- 
tries in mixed final states. In principle, these measure- 
ments allow us to determine all the relative sizes and 



phases of the amplitudes A K and A K , which are related 
to <f>2 through an isospin relation 0, @| by 



D +2i 



A+ + A-+ 2A° 
A+ + A-+ 2A° 



(13) 



Consequently, in the limit of high statistics, we can con- 
strain 4>i without discrete ambiguities. 



To reconstruct candidate B° 



decays, we 



combine pairs of oppositely charged tracks with 7r° can- 
didates. The selection criteria for charged tracks are the 
same as in the previous B° — > p ±f K T analysis Can- 
didate 7r°'s are reconstructed from 7 pairs having M 77 
in the range 0.1178-0.1502 GeV/c 2 , corresponding to ±3 
standard deviations (a) in M„o resolution, and momenta 
greater than 0.1 GeV/c in the laboratory frame. We 
require E 1 > 0.05 (0.1) GeV in the barrel (endcap) of 
the electromagnetic calorimeter Q, which subtends 32°- 
129° (17°-32° and 129°-150°) with respect to the beam 
axis. Candidate B mesons are reconstructed using two 
variables calculated in the center-of-mass frame: the B 
invariant mass calculated using the beam energy in place 
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of the reconstructed energy (Mb c ), and the energy dif- 
ference between the B candidate and the beam energy 
(AE). We define a signal region -0.1 GcV < AE < 
0.08 GcV and 5.27GcV/c 2 < M hc , and a large fitting 
region \ AE\ < 0.2 GeV and 5.2 GeV/c 2 < M bc . 

The procedure used to measure At and to determine 
the flavor of the decaying B° meson, q t ag, and its qual- 
ity, I, are described elsewhere [lij]. The dominant back- 
ground is e + e~ — > qq (q = u,d,s,c) continuum events. 
To distinguish these jet-like events from the spherical B 
decay signal events, we combine modified Fox- Wolfram 
moments [l5| and the B flight angle with respect to the 
beam direction into a signal (background) likelihood vari- 
able £ s i g (bkg) and impose requirements on the likelihood 
ratio 1Z = £ S i g /(£ S i g + £bkg)- These requirements de- 
pend on the quality of flavor tagging. When more than 
one candidate in the same event is found in the large fit- 
ting region, we select the best candidate using likelihood 
based on M 77 and 1Z. After the best candidate selection, 
we apply a Dalitz plot cut: candidates are required to 
satisfy 0.55 GeV/c 2 < ^/s^J < 1.0 (0.95) GeV/c 2 for at 
least one of s+, s_, or sq. In the fits below, we use square 
Dalitz plot variables (m', 9') for convenience, performing 
a parameter transformation on (s+, s_) [16j . 

Figure [1] shows the Mb c (AE) distribution for the re- 



constructed B° 



7T candidates within the AE 



(-Wbc) signal region. The signal yield is determined 
from an unbinned four-dimensional extended-maximum- 
likclihood fit to the AE-M bc and Dalitz plot distribution 
in the large fitting region, where the Dalitz plot distri- 
bution is fitted only for events inside the AE-M bc signal 
region. The fit function includes signal; incorrectly re- 
constructed signal, which we call self-cross- feed (SCF); 
continuum: and BB background components. The prob- 
ability density function (PDF) for each component is the 
same as that used for the TDPA described below, but 
integrated over At and summed over qt&g- The fit yields 
971±42 B° — > 7T + 7r~7r° events in the signal region, where 
the errors are statistical only. 

Using the same data sample as described above but 



>300 
^250 




[(b) 


— signal 




SCF 




BB 




A ^qq 




W \ 


1 






WWa.... 



5.22 5.24 5.26 5.28 5.3 

M bc (GeV/c 2 ) 



0.2 



-0.1 0.1 

AE (GeV) 



0.2 



FIG. 1: (a) A/bc and (b) AE distributions within the AE 
and Mb c signal regions. Solid, dot-dashed, dotted, and dashed 
hatched histograms correspond to correctly reconstructed sig- 
nal, SCF, BB, and continuum PDF's, respectively. 



performing a time-integrated Dalitz plot fit with a wider 
Dalitz plot acceptance, 0.0 (0.55)GeV/c 2 < ^/s±(o) < 
1.5GeV/c 2 , we determine the p lincshape parameters (3 
and 7. We use the results obtained for the TDPA be- 
low. We also put upper limits on the possible deviations 
of ( /3 K , { 7^) from the nominal (/3, 7), which we use to 
estimate systematic errors. 

To determine the 26 coefficients, we define the follow- 
ing event-by-event PDF: 

P= J2 fxVx(AE, M bc , m', 9', At, g tag , , 

X=s\g,qq,B~B 

where V s i E , V q -q, and V b -q are the PDF's of signal includ- 
ing SCF, continuum, and BB components, respectively, 
and / s ig, f q q, and / B -g are the corresponding fractions 



that satisfy / sig + f qq + f 



'BB 



1. Here, Psie and V 



BB 



arc modeled based on Monte Carlo (MC), though a small 
correction is applied to 'Psig to take account of the differ- 
ence between data and MC, while V qq is modeled using 
data. The signal PDF, 'Psig, is the sum of a correctly 
reconstructed PDF (Vt rne ) and an SCF PDF, where 

^true *^true^true 

(AE,M hc )e(m',6';l) 
rue (m! ,9' ,At,q taE ;l) ■ 

Here ^ ruc , V tmc (AE , M hc ) , e(m',9';l), and 
^trueCwi'j 8'> At, (frag! are event fractions in each 
category of tagging quality I, a AE-M bc PDF, a Dalitz 
plot dependent efficiency, and a Dalitz-At PDF for the 
correctly reconstructed signal component, respectively. 
The Dalitz-Ai PDF corresponds to the right-hand side 
of Eq. (fl]) with the following modifications: (i) it is 
convolved with the At resolution function [T3|; (ii) it is 
multiplied by the determinant of the Jacobian for the 
transformation (s+,s_) i— > (m',9 1 ); and (hi) the wrong 
tag fractions, wi, and the difference between B° and B° 
decays, Awi, are taken into account. A more detailed 
description of the PDF can be found elsewhere [l8j |. 

An unbinncd-maximum-likclihood fit to the 2824 
events in the signal region yields the results listed in Ta- 
ble |TJ With a toy MC study, we find that the errors 
estimated by the likelihood function do not give correct 
68.3% confidence level (C.L.) coverage for the interfer- 
ing parameters. In the table, we multiply the error es- 
timates from the likelihood function by a factor of 1.17, 
which is calculated from the MC study, to obtain errors 
with correct coverage. We find that Uq is 4.8 a above 
zero, corresponding to clear evidence for the presence of 
the decay B° — * p°ir° in agreement with our previous 
measurement [nj. Figure [H shows the mass and helicity 
distributions, and the background-subtracted At asym- 
metry plot for each pix enhanced region. We define the 
asymmetry in each At bin by (N+ - N-)/(N+ + N-), 
where corresponds to the background-subtracted 

number of events with <7 ta g = +1 ( — 1)- The p~n + en- 
hanced region shows a significant cosine-like asymmetry 



5 



TABLE I: Results of the time-dependent Dalitz fit (left three columns), 
parameters (rightmost column), whose definitions can be found elsewhere | 
systematic, respectively. The correlation coefficient between At" and ApJ~ 



and the associated quasi-two-body CP violation 
The first and second errors are statistical and 
(Apo^o and Spt^o) is +0.47 (-0.08). 





+ 1 (fixed) 


u+ 


+0.23 + 0.15 + 0.07 


1+ 


-0.01 + 0.11 + 0.04 


U+ 


+f.27±0.f3±0.09 


uz 


-0.62 + 0.16 + 0.08 


7_ 


+0.09 + 0.10 + 0.04 




+0.29 ±0.05 ±0.04 




+0.15 + 0.11+0.08 


la 


+0.02 + 0.09 + 0.05 




+0.49 + 0.86 + 0.52 




-1.18 + 1.61 + 0.72 


rRo 
■* + - 


+1.21 + 2.59 + 0.98 


rr+,Re 
U +0 


+0.29 + 0.50 + 0.35 


rr-,Re 


-2.37 + 1.36 + 0.60 


r Re 

y +o 


+1.15 + 2.26 + 0.92 


j j+, Re 
U -0 


+0.25 + 0.60 + 0.33 


77-Re 


-0.53 + 1.44 + 0.65 


rRe 

J--0 


-0.92 + 1.34 + 0.80 


v +,lm 


+1.18 ± 0.86 ±0.34 


UZ'l m 


-2.32 + 1.74 + 0.91 


rim 

■* + - 


-1.93 + 2.39 + 0.89 




-0.57 ±0.35 ±0.51 




-0.41 ± 1.00 + 0.47 


rim 
^+0 


-0.40 + 1.86 + 0.85 




-1.34 ±0.60 ±0.47 


r/-' Im 

U -0 


-0.02 ± 1.31 ±0.83 


rim 

■*-0 


-2.03 + 1.62 + 0.81 



C 

AC 
S 

AS 



-0.12 + 0.05 + 0.04 
-0.13 + 0.09 + 0.05 
+0.36 + 0.10 + 0.05 
+0.06 + 0.13 + 0.05 
-0.08 + 0.13 + 0.05 



+0.21 + 0.08 + 0.04 
+0.08 + 0.16 + 0.11 



A p a n o 



-0.49 ± 0.36 ± 0.28 
+0.17 + 0.57 + 0.35 



arising from a non-zero value of UZ . Note that this is not 
a CP-violating effect, since p~n + is not a CP-eigenstate. 
No sine-like asymmetry is observed in any of the regions 
(g)-(i)- 

The non-interfering parameters can be interpreted as 
the quasi-two-body parameters of the process B° — > 
p^ir^ 1 , whose definitions can be found elsewhere [H}, and 
the CP violation parameters of the process B° — > p°ir°: 
ApO^o = —UZ/Uq and S p o v o = 2Iq/Uq. These are also 
listed in the Table [0 

There are several sources of systematic uncertainty. To 
determine their magnitudes, we vary each possible con- 
tribution to the systematic error by its uncertainty in the 
data fit or in the MC, and take the resultant deviations in 
the fitted parameters as errors. We add each contribution 
in quadrature to obtain the total systematic uncertainty. 
The largest contribution for the interfering parameters 
comes from radial excitations. We take account of pos- 




0.6 0.8 1 -1 1 -8 8 

m K (GeV/c ) cos9 K At (ps) 



FIG. 2: Mass (a)-(c) and helicity (d)-(f) distributions, and 
background subtracted At asymmetry plots in the good tag- 
ging quality region I > 3 [14| (g)-(i), corresponding to the 
p+n- [(a),(d),(g)], P -7T+ [(b),(e),(h)], and p°n° [(c),(f),(i)] 
enhanced regions. The notations for histograms (a)-(f) are 
the same as Fig. Q] 



sible deviations of ( Y K ) from the (/3, 7) values, and 
uncertainties of (3, 7, and the mass and width of each res- 
onance. Large contributions to the systematic errors for 
the non-interfering parameters come from potential back- 
grounds such as B° -»■ /o(980)7r ,/ (600)7r°,a;7r , and 
non-resonant 7r + 7r _ 7r°, which we neglect in our nominal 
fit. We perform fits to toy MC including these back- 
grounds with the branching fractions at their 68.3% C.L. 
upper limits, which we obtain from our data or world 
averages (IH l2fl | ; the largest variations are taken as sys- 
tematic errors. Comparable contributions also come from 
vertex reconstruction, background PDF's, and tag-side 
interference [2l[ ; more detail can be found elsewhere [l8| . 

We constrain </> 2 from the 26 parameters measured in 
our analysis following the formalism of Ref. [H| and the 
statistical treatment using toy MC described in Ref. [13] • 
The resulting 1 — C.L. function is shown in Fig. [3] as a 
dotted curve. To incorporate all available knowledge, we 
combine our measurement with results on the branch- 
ing fractions for B° — ► p ± 7r T and B + —+ _P + tt°, P°tt + , 
and flavor asymmetries of the latter two [20(. Assuming 
isospin (pentagon) relations @, 0] and following the same 
procedure as above, we perform a full Dalitz and pen- 
tagon combined analysis, the result of which is shown in 
Fig. [3] as the solid curve. We obtain 68° < 4> 2 < 95° as 
the 68.3% confidence interval for the solution consistent 











\ ^ C.L.=68.3% 



30 60 90 120 150 180 
<Sf 2 (degrees) 

FIG. 3: 1 — C.L. vs. (j>2- Dotted and solid curves correspond 
to the result from the TDPA only and that from the TDPA 
and an isospin (pentagon) combined analysis, respectively. 
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with the SM expectation. A large SM-disfavored region 
(0° < 2 < 5°, 23° < 2 < 34°, and 109° < 2 < 180°) 
also remains. In principle, with more data we may be 
able to remove the additional 02 solutions. 

In summary, using 4141b" 1 of data we have performed 
a full Dalitz plot analysis of the B° — > ir + ir~ir° decay 
mode, where the observables include the first measure- 
ment of SpO^o. A full time-dependent Dalitz plot analysis 
with the pentagon isospin relation is performed for the 
first time and a constraint on the angle 02 is obtained. 
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